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Through the use of a lattice U(l) Ward-Takahashi identity, one can find a precise 
definition of flux and electric four-current that does not rely on the continuum 
limit. The magnetic four-current defined for example by the DeGrand-Toussaint 
construction introduces order a 2 errors in the field distributions. We advocate 
using a single definition of flux in order to be consistent with both the electric 
and magnetic Maxwell's equations at any lattice spacing. In a U(l) theory the 
monopoles are slightly smeared by this choice, i.e. are no longer associated with a 
single lattice cube. In Abelian projected SU(2) the consistent definition suggests 
further modifications. For simulations in the scaling window, we do not foresee 
large changes in the standard analysis of the dual Abrikosov vortex in the maximal 
Abelian gauge because the order a 2 corrections have small fluctuations and tend 
to cancel out. However in other gauges, the consequences of our definitions could 
lead to large effects which may help in understanding the choice of gauge. We also 
examine the effect of truncating all monopoles except for the dominant cluster on 
the profile of the dual Abrikosov vortex. 



1. Introduction 

Dual superconductivity has long been suggested as a possible mechanism 
for quark confinement signaled by a spontaneously broken U(l) gauge sym- 
metry and manifested by a dual Abrikosov vortex between quark and anti 
quark. This picture was verified some time ago in Abelian projected SU(2) 
lattice gauge theory in the maximal Abelian gauge 1 . More recently further 
studies have elaborated this picture 2,3,4 . 
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As in all lattice calculations, there is freedom in choosing lattice opera- 
tors, requiring only that they agree with the continuum definition to lowest 
order in the lattice spacing a. However we have the opportunity in these 
studies to be more precise by incorporating the lattice Ward-Takahashi 
identity derived from the residual U(l) gauge symmetry 5 . This gives an 
Ehrenfest relation for the expectation value of the fields and currents giving 
the electric Maxwell equations exactly at finite lattice spacing. Interest- 
ingly, this defines a unique lattice expression for the field strength or flux 
and the electric and magnetic currents. 

In the present work, we examine the impact of this on the study of the 
dual Abrikosov vortex. To be consistent, the magnetic Maxwell equation 
must use the same definition as the electric one. However the standard 
DeGrand-Toussaint 6 [DT] definition of the magnetic current is based on 
a different definition of flux, resulting in inconsistencies in the magnetic 
Maxwell equation. We argue here that one should alter the DT construc- 
tion, using a single definition of flux throughout. A consequence is that the 
magnetic current no longer contains discrete monopoles but rather a more 
general magnetic charge distribution. In effect, the monopoles are smeared 
in our picture. 

This consistency question is only relevant at finite lattice spacing and all 
these concerns go away in the continuum limit. Nevertheless it is desirable 
to have a consistent treatment of total flux in the vortex determined from 
the electric Maxwell equation, and the profile shape determined by the 
magnetic Maxwell eqution at fixed lattice spacing. Further we note that 
the finite lattice spacing effects are significant for the values of (3 that we 
often use for calculations. 

We also report on the effect of truncating DT monopole loops, keeping 
only the one large connected cluster. This truncation is expected to have no 
effect on the confinement signal 9 ' 10 . This should manifest itself here in that 
the tail of the profile of the magnetic current of the vortex are unaffected 
by the truncation. We find this to be the case. This procedure requires 
that the magnetic current consists of discrete monopoles. Truncation is 
not well defined in the above smeared monopole picture. Hence we take 
the conventional view in presenting these results that the continuum limit 
is required to obtain the consistency described above. 
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2. Three definitions of flux 

Let us consider three definitions of field strength, all agreeing to lowest 
order in a. The first definition was used by DT to define monopoles: 

F$ = e^m) - M (m + i/) - e v {m) + e v {m + u) - 
= 6n V (m) - 2nn^, v , 

where M (m) refers to the U(l) link angle in the domain — n < 8^ < +ir. 
The integers n^ v are determined by requiring that — 7r < F$ < +tt. That 
is F$ is a periodic function of 6^ v with period 2tt. Here quantities with 
mean those which appear in the lattice numerical calculation without 
appending factors of e and a. 

The second and third definitions gives the exact electric Maxwell equa- 
tion for lattice averages 

A-(F$) w = (j;U) w , i = 2,3, (1) 

where 

{ "' )w= (e»") ' 

for the cases of U (1) gauge theory and SU (2) gauge theory respectively. We 
give the U(l) derivation in detail since it is straight forward and contains 
the essential points of the argument. There are significant complications in 
the SU (2) case and hence we just sketch that derivation. 

2.1. Flux in the U(l) gauge theory 

Consider 

Z w (e„(m)) = J [d9] sin9 w exp (J3S, ) 

S = ^2 cos6» M „(n), 0=~2- 

The subscript of Zn/(e /J (m)) refers to the incorporation of the source into 
the partition function and the argument is a variable defined as the shift 
of one particular link, 9^(m) — > M (m) + e^m). This translation can be 
transformed away since the measure is invariant under such an operation. 
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Therefore 

= 5Z W = [ [d9] sin 9 W exp(/35) - f [d9] sin 9 W exp(pS) 

= e M J [dO] (d^ (m) cos 9 W + sin 9 W exp (f3S) , (2) 

where 5^(m) = ±1 if m labels a ± directed link and = otherwise. This 
is the static current generated by the Wilson loop. 

5,{m) = j; (3) 

Next evaluate the derivative of S. Isolate the six plaquettes affected by 
the shift 

S = [cos{6^(m) + e v (m + fi) -9 ll {m + v)- 6 v {m)} 

+ cos{6' At (m) - 9 v (m + fi—v) — 9^(m - v) + 6 v (m - v)}] H 

dS 

Y [- sin {(9 M (m) + 6 u {m + -O^m + v) - v {m)} 



de^m) 

- sin {6fj,(m) - 9 v {m + y, - v) - 9^(m — v) + 9 v (m - v)}] , 
dS 

^_ y = -A;sin^(m). (4) 

Using Eqn.(l) and Eqn.(4) we can see that Eqn.(2) is the form of Maxwell 
equations for averages. 

1 



;K{F$)w = J e Jm), (5) 



e 

where 

(sm9 w ---) 
(cos 9 W ) 

Since the charged line in a Wilson loop is closed the electric current is 
conserved. The local statement of conservation is 

q = a;a-(f$) w = a-j;. 

It is straightforward to verify on the lattice that the LHS of Eqn.(5) gives 
— 1,0,1 depending on its position with respect to the Wilson loop. An 
alternative definition such as F$) need not vanish off the Wilson loop nor 
give ±1 on the Wilson loop and hence would introduce an error. In a U(l) 
theory there is no dynamical charge density, all charge resides on the Wilson 
loop. 



February 1, 2008 4:20 WSPC/Trim Size: 9in x 6in for Proceedings 



c2003procvl 



5 



2.2. U(l) flux in the SU(2) theory in the maximal Abelian 
gauge 

We restrict our attention to the maximal Abelian gauge defined as a local 
maximum of 

R = J2^{^Mn)a 3 Ul(n)}, 

over the set of gauge transformations {g(m) — e la i( m ) a i\ jj — > JJ9_ Tak- 
ing U to be the stationary value, the stationary condition is given by 

dRp«] 



F jn [U] 



= 0. 

a=0 



Jni ~ i da 3 (n) 
The second derivatives entering in the Jacobian are given by 



M- - (U) d2Rm 



a=0 



daj(n)dai(m) 
The partition function is 

Z 9 w f -(el(m)) = f [dU] ^Tr[ia s U w (n)] e X p((3S)l[S(F jn [U}) A FP ,(6) 

jn 

where the Faddeev-Popov Jacobian is 

A FP = det\M jn;im (U)\. 

An infinitesimal shift in this partition function has the added complication 
that it violates the gauge condition. This can be corrected by an infinitesi- 
mal accompanying gauge transformation. Thus the shift in one link affects 
all links. However experience has shown that the effect drops off rapidly 
with distance from the shifted link. 

The derivative of the partition function Eqn.(6)with respect to e^(m) 
gives 



= J [dU] h»(rn)^Tr[U w (n)]\exp(f3S)l[6(F jn [U])A FP 

jn 

+ J [dU] |^lTr[ i a 3 ^(n)]^-|^|exp(/35)n^(^[t / ])AFP 
+ J [dU] l -Tr^U w {n)]e^{pS) W[5{F jn [U\)* F p \ . (7) 

M [ jn J 

The third integral contains terms in the Ward Takahasni identity coming 
from the gauge fixing including ghost contributions. 



February 1, 2008 4:20 WSPC/Trim Size: 9in x 6in for Proceedings 



c2003procvl 



6 



where 



We can cast this into the form of the electric Maxwell's equations for 
averages as in the case of the U(l) theory. However there are now more 
terms in the current. Consider the standard U(l) parametrization of an 
SU(2) element: 

( C^(n)e ie ^ ^(n)e^( n )-^( n )>\ 
'„ S^\(e i0 - 

C^n) = cos^(n), 

S^(n) = sin^(n). (8) 

In the Abelian projection factored form, the righthand factor contains the 
1/(1) photon, parameterized by 0. The lefthand factor contains the charged 
coset matter fields, parameterized by <p and 7. The transformation proper- 
ties are well known and reviewed in DiCecio et.al. 5 . 

We consider an alternative separation into diagonal and off-diagonal 
parts which is needed in defining the flux. 

U ^ n >-{ C ll e-^) + \-S lx e- i ^-°^ )' 
= £> M (n) + M (n). 

The off-diagonal part is the charged matter field 4?^ = S^e - ^ 7 ^ -9 ^. and 
diagonal part includes the photon, e 10 *, but also a neutral remnant of the 
matter field a/1 — |$ M | 2 which — ► 1 in the limit a — > 0. 

To cast Eqn.(7) into the form of a current conservation law, we first 
consider the terms to zeroth order in M . First the Wilson loop. Isolating 
the diagonal contributions gives 

Uw = Dw + Uw> 
where D\y is the product of the diagonal parts 

^Tr[D w (n)] = JJ(7 M (n) • cos8 w , 

w 

^Tr{icr 3 D w (n)} = - Y[C^n) • sin W . 



w 
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Wc adhere to the standard choice of an Abelian Wilson loop in which we 
drop any contributions due to off diagonal elements and further take 
the factors C^n) = 1 giving 



2 

Second, consider the action. Write 



irr[C/Abelian (n)] = cos ^ ; 
Tr[ia 3 U& beli&n (n)] = -sm9 w . 



S = E ^[D„An)} + S, 



2 

where S contains terms involving 0^{n). 
1 



d(S-S) _ x - 
deUm) ~ 2f 



2 Tr {-ia 3 (m)Dp(m)D v (m + /i)Dj,(m + v)Dl(m)} 



+ ^Tr {-ia 3 (m)D^(m)Dt(m + n~ f)£>t(m - v)D v {m - v)} 
Since all matrices are diagonal we can simplify: 

^T7~T = S A ^ [C M (m)a(m + /i)C„(m + i/)a(m) sin^(m)]. 

The quantity in square brackets is antisymmetric in /x^ and we identify this 
as proportional to the field tensor. 

F$ = Cf, (m)C v (m + (j,)Cn(m + is)C v (m) sin M1/ (m) 

Returning to the identity, Eqn.(7), and using the notation 

^ [d£/]{. • • } exp (/3S) [] <5(^„[C/])A FP = (• • -) g . f . 

we obtain 



= <S M (m)(cos0 w > 9 ./. ~ P(^0 W A-F^) g . f . - p(sme w g ^^ ) g . f . 



+ gauge fixing terms + ghost terms . 
Rearranging terms as in the U(l) case, we get 

dS 



(sme w AzF$) g . f . = s (m) _ ^ sin ^a^m))^ + g.f. fe ghosts 
(cos0w) s ./. m (cos6 w ) g .f. (cos6w}g.f. , 
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K(F$)w, g .f. = 4 e) Abelian Wilson loop + h {e) matter fields 



+ 



J M (e) gauge fixing + j^(e) ghosts ^ 

The 'Abelian Wilson loop' term is analogous to the above U(l) case. The 
'charged matter field' term arises from the off-diagonal elements of links 
in the action expression and would contribute without gauge fixing. The 
'gauge fixing' term arises from the corrective gauge tranformation that ac- 
companies the shift of a link. The 'ghost' term arises from the shift and 
accompanying gauge transformation on the Faddccv-Popov determinant. 
See DiCecio et.al. 5 for a complete derivaton of all terms. 

It is important here to note that we are not interested in distinguishing 
the various contributions to the current in the present work. We are only 
interested in the total current and that can be obtained from the LHS of 
Eqn.(9). 

3. Consistency with the magnetic Maxwell equation 

Having defined a unique flux F$ through the electric Maxwell equation, 
the magnetic Maxwell equation is 

-f a+fW — f m ( i ) _ o Q 

2 ^iivpa^u 1 pa ~ J p l — A,o. 

However the standard DT definition of current is 



rm(l) 1 A+F 



(i) 

pa ■ 



and hence if we use the conventional to define the monopole current, 
and or F^ respectively for U(l) and SU{2) theories to get an ex- 
act expression for flux in the confining string, then the magnetic Maxwell 
equation is violated. 

The solution is to relax the requirement that we use the DT monopole 
definition and use F^ or F^> instead when defining monopoles. A simple 
configuration for the U(l) case (F^) will illustrate the effect. Consider a 
single DT monopole with equal flux out of the six faces of the cube. Then 
the ratio of the F^ 2 > flux out of this cube compared to the DT F^ 1 ' flux 
gives 

6 7^ « 0.83. 
6(2tt/6) 

Since the current is conserved, the balance is made up by magnetic charge 
in the neighboring cells. On a large surface the total flux is the same for 
the two definitions. 
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In Fig.(l) wc plot the F$ as a function of 9^ v , giving a "sawtooth" 
shape. Monopoles occur as a consequence of F^ 1 ' crossing the sawtooth 
edge, giving a mismatch of 2ir in the flux out of a cube. The sine function 
has no such discreteness and so the notion of discrete Dirac strings and 
Dirac monopoles is modified. However as one approaches the continuum 
limit, the action will drive the plaquette to zero and then the regions where 
the sawtooth differs from the sine are suppressed. Hence we expect both 
forms to give the standard Dirac picture in the continuum limit. 




Figure 1. F$ (sawtooth) and F$ (sine) as a function of the plaquette 
angle 9^ 

The electric Maxwell equation determines the total flux in the confining 
string and the magnetic Maxwell equation determines the transverse shape. 
Further the latter enters directly in the determination of the London pen- 
etration length, Ad- To see this let us consider the classical Higgs theory 
which we use to model the simulation data. The dual field is given by 

G^{m) = A-ei d \m) - A^m), 

where 6jf'(m) is a dual link variable. Let us choose to break the gauge 
symmetry spontaneously through a constrained Higgs field. 

$(m) = vp(m)e lx( - m \ p(m) = 1. 

Under these conditions the magnetic current simplifies to 

J™(™) = 2e r 2 „w 2 sin{^(m) + x M (m + /z) -^(m)}. 
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where e m is the magnetic coupling. The magnetic Maxwell equation is 

A+G„„ = J?. (10) 

For small 0^ it is easy to see that there is a London relation of the 
form 

d ^ {m) = 2^y ( A m ^(H - A" J™(m)) . (11) 

Taking the confining string along the 3 axis and choosing \i — 1 and v = 2 
we see that the profile of the third component of curl of the magnetic current 
must match the third component of the electric flux profile. This assumes 
an infinite Higgs mass Mh- With a finite mass there is a transition region 
of size ~ 1/Mjj in the core of the vortex but the above London relation 
must hold sufficiently far outside the core. 

Combining Eqns.(lO) and (11) we get the relation 

(1 - A^A+A;) (E 3 (m)) w = 0, A- 1 = e m vV2 

The corresponding equations in the simulation of the SU{2) theory must 
also be satisfied in order to arrive at this correct expression for A<j and 
hence the importance of our definitions. 

4. Numerical Results 

We generated 208 configurations on a 32 4 lattice at (3 — 2.5115. The maxi- 
mal Abelian gauge fixing used over- relaxation. Fig. (2) shows the profile of 
the electric flux corresponding to and F^ 3 \ Fig. (3) shows the profile 
of the theta component of the magnetic current corresponding to F^ 
and F< 3 ). 

In summary we showed that consistency requires one use the same 
definition of flux throughout. If, for example, one uses F^ definition of 
electric field (bottom graph in Fig. (2)) in order to account correctly for 
the total flux but then uses the DT definition of current (top graph in 
Fig. (3)) we would then incur errors of ~ 40%. For simulations in the 
scaling window, as long as we use the same definition consistently we do 
not foresee large changes in the standard analysis of the dual Abrikosov 
vortex in the maximal Abelian gauge because the order o 2 corrections have 
small fluctuations and tend to cancel out 8 . However in other gauges, the 
consequences of our definitions could lead to large effects which may help 
in understanding the choice of gauge. 
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Figure 2. Profile of the electric field (highest to lowest) E$ (circles) and 
£3 (diamonds) on the mid-plane between q and q separated by 13a. 




Figure 3. Profile of the theta component of the magnetic current on the 
mid-plane between q and q separated by f 3a based on (highest to lowest) 



F^J (triangles), F\fJ (diamonds) and F^f (inverted triangles). 

Finally we report on the effect of truncating DT monopole loops, keeping 
only the one large connected cluster 9 ' 10 . As mentioned in the Introduction 
this truncation is expected to have no effect on the confinement signal. This 
should manifest itself here in that the tail of the profile of the electric field 



?(2) 



(3) 
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and magnetic current of the vortex are unaffected by the truncation. This 
is born out as expected. Fig. (4) shows that radial profile of the magnetic 
current is indistiguishable except for a small deviation in the core of the 
vortex. 

0.1 i— 




4 8 12 

radial dist. 



Figure 4. Profile of the theta component of the magnetic current on the 
mid-plane between q and q separated by 13a based on (highest to lowest) 
FjxJ (circles), F$ truncated (diamonds). 
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